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ABSTRACT 

We analyze the low-redshift (z«0.2) circumgalactic medium by comparing absorption¬ 
line data from the COS-Halos Survey to absorption around a matched galaxy sample 
from two cosmological hydrodynamic simulations. The models include different pre¬ 
scriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), 
and constant winds (cw). We extract for comparison direct observables including 
equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are 
generally in good agreement with these observations for HI and certain low ionization 
metal lines, but show poorer agreement with higher ionization metal lines including 
Sim and Ovi that are well-observed by COS-Halos. These discrepancies suggest that 
both wind models predict too much cool, metal-enriched gas and not enough hot gas, 
and/or that the metals are not sufficiently well-mixed. This may reflect our model as¬ 
sumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a 
heuristic prescription to quench massive galaxies by super-heating its ISM gas, which 
we show yields sufficient low ionisation absorption to be broadly consistent with obser¬ 
vations, but also substantial O VI absorption that is inconsistent with data, suggesting 
that gas around quenched galaxies in the real Universe does not cool. At impact pa¬ 
rameters of < 50 kpc, recycling winds dominate the absorption of low ions and even 
Hi, while Ovi almost always arises from metals ejected longer than 1 Gyr ago. The 
similarity between the wind models is surprising, since we show that they differ sub¬ 
stantially in their predicted amount and phase distribution of halo gas. We show that 
this similarity owes mainly to our comparison here at fixed stellar mass rather than at 
fixed halo mass in our previous works, which suggests that CGM properties are more 
closely tied to the stellar mass of galaxies rather than halo mass. 


1 INTRODUCTION 

Studying the circumgalactic medium (CGM) is an emerg¬ 
ing frontier of galaxy evolution that connects galaxies to 
the intergalactic medium. In this tenuous multi-phase gas 
that surrounds galaxies lie clues to the processes that 
drive gala xy evolution such as gaseous inflows and out¬ 


flows (e.g. Kereset_al, 1 120051 : iDekel et al.l 120091 ; iDave et all 
120121 ; iLillv et al.ll2013h . Recent observations with Hubble’s 
Cosmic Origin Spectrograph (COS) have characterised the 
CGM in unprecedented detail. They have shown that there 


is probably more oxygen in t he halos around star-f orming 
galaxies than within them dTumlinson et all 1201 il l . that 
the CGM may contain th e missing baryons and missing 
meta l s within galaxy ha los llPeeples et al .1120141 : IShull et al.l 
I20l4 IWerk et all 120141 1 . and that star-f orming and pas¬ 
sive g alaxies have different Ovi signatures jTumhnsone£_al. 
2011) but similar cool gas cont ent ( Thom et all 120121 ; 


Wer^_et_aL Il20l4). R ecent work bv IStocke et all (120131) and 
Tripp &; Sond ( 20121 ) has also used COS to probe the CGM 


out to several hundred kiloparsecs, and observations by 
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IProchaska et al.l d201 ll'l and lRudie et al.l d2013h have estab¬ 
lished a characteristic extent of the metal-enriched CGM 
to be roughly ~ 300 proper kpc, regardless of redshift. 
In addition to these observations that characterise the 
CGM, there is a wealth of obser vational data d etecting out¬ 
flows at a variety of redshifts 
2005| Ttemonti etaL 2007; 


hifts llM artinl 120051: iRupke et al.l 

; Weiner 200f^ Rubinet^L 20121: 


Pettini et al.l 200 ll : Steidel l200ll ; Veilleux et al.l 2005ll . indi¬ 
cating that most if not all galaxies drive winds at some point 
in their evolution, presumably enriching the CGM and IGM 
with mass and metals in the process. 

The physical mechanisms that drive these outflows are 
not well characterised. While observations can now probe 
individual sight lines around individual galaxies, there is 
not yet a robust understanding of the metallicity, covering 
factor, and thermodynamic phase of the outflowing mate¬ 
rial. The specifics of how these quantities vary with the 
galaxy’s mass, star formation rate, environment, and red- 
shift are poorly characterised. The fate of the ejected mate¬ 
rial after it leaves the galaxy is even less well constrained. 
Simulations suggest the presence of “halo fountains”, where 
baryons are ejected from th e galactic di sk only to rain back 
down again at later times lOppenheimer fe Davdl2008ll . The 
contribution of this so-called recycled accretion relative to 
pristine inflo ws from the IGM depen ds at minimum on 
gala xy mass (lOp penhei mer et al.ll2010tl and impact param¬ 
eter llFord et a! 20l4l . and plausibly also on star formation 
rate (SF R), SFR surface density, a nd environment. Simula¬ 
tions by lOppenheimer et ali (j2012l l suggested that the IGM 
is enriched “outside-in”, where winds propagate further from 
galaxies at earlier epochs and metals re-collect towards halos 
at late times, but this conclusion may depend on the poorly- 
understood mechanisms by which winds are launched. 

Hydrodynamic simulations that directly model inflows 
and outflows provide a crucial complement to emerging 
CGM observations by plausibly elucidating how the physical 
conditions and dynamical st ate of the CGM ga s can be man¬ 
ifest in absorption line data. iFord et aD (|2013h showed that 
metal ions with low ionisation potential (e.g., Mg II) trace 
dense gas close to galaxies, while metals with high ionisa¬ 
tion potentials (e.g., Ovi) trace more diffuse gas. This work 
also showed how different prescriptions for outflows can re¬ 
sult in significant differences in th e temperature , density, 
and metallicity of the CGM gas. IFord et all d2014|j followed 
and showed, in detail, how the baryon cycle operates in a 
plausible wind model. At low redshift, inflowing material is 
dominated by recycled gas that was previously ejected in 
an outflow; hence inflows are enriched and can be traced 
by metal lines. In general, low metal ions like Mgn trace 
this recycled accretion, while Ovi traces hotter, less dense 
material built up by prior epochs of outflows. 

While these simulations have provided some general 
guidelines for interpreting CGM observations, they have yet 
to be tested in detail against existing CGM data. The best 
existing absorption line data set that targets a range of 
impact parameters close to galaxie s in the CGM is cur¬ 
rently from the COS-Halos survey (iTumlinson et all l201ll . 
l2013l : IWerk et ai1l2012l . l2013ll . This survey targets 39 quasars 
around 44 0.1L* — 3 L* galaxies with a carefully-selected 
sampling of colour and stellar mass, and constitutes the 
first galaxy-selected absorption line survey of low-redshift 
CGM gas. These galaxies have stellar masses ranging from 


10 9 , 5 ~ii- 5 ) with spectroscopic redshifts z v hot « 0.14 — 0.34, 
and are located 15 — 160 kpc (projected) from a background 
UV-bright quasar. While COS-Halos covers a wide range of 
ions, for this work we focus on Hll216, Silll260, Clll335, 
SiIII1206, Siivl394, Ciii 977, and Ovi, since they span a 
representative range of ionisation energies (see Table 0. 
Typical COS-Halos detection thresholds for the various ions 
are given in Table [l] and are generally around rest equiv¬ 
alent widths of O.lA. The COS-Halos sight lines are sup¬ 
plemented by Keck/HIRES echelle spectra for 35 quasars to 
include coverage for Mgn (2796, 2803A). These observations 
represent the largest sample of CGM absorption line data at 
these impact parameters, and hence provide a stringent test 
for models of CGM gas and their relationship to galaxies. 

In this paper we build on our previous works us¬ 
ing hydrodynamic simulations to (i) determine if our hy¬ 
drodynamic simulations quantitatively reproduce the COS- 
Halos absorption line observations; (ii) test whether such 
observations can distinguish between plausible but sub¬ 
stantively different models for galactic outflows; and (iii) 
constrain the physical and dynamical state of CGM gas 
around ~ L* galaxies as probed by COS-Halos data. 
As observati ons improve with o ngoing surveys such as 
COS-Dwarfs (iBo rdoloi et al.11201 41), the COS-GTO survey 


llPanforth et al.1 l 201 p|^ 20 U, 2014; IStocke et al.l l2013l l , and 
other efforts (e.g. lHeckma n l2013l ) to study the CGM of se¬ 
lected galaxy samples, this work aims to serve as a template 
for how CGM absorption line observations can constrain the 
physical processes driving the baryon cycle. 

Our paper is organised as follows: in §2 we describe our 
simulations, in §3 we compare to direct observables, in §4 
we describe variations with wind model including variations 
in baryon cycling, and in §5 we summarise our conclusions. 


2 SIMULATIONS AND ARTIFICIAL SPECTRA 
2.1 The code and input physics 

For this paper we use our modified version of the N- 
body+entropy-conserving s mooth particl e hydrodynamic 
(E C-SPH) code Gadget - 2 (tSpringell 120051 1. fully described 
in lOppenheimer fc Pavel ([2008). Our main simulation is 


identi cal to that used in Pave et ~aT1 (l2013l l: IFord et aD 
d2014l') . and further details can be found there. Briefly, 
we assume a WMAP-9 concordant ACPM cosmology 
(iHinshaw et al.l 120091 1 : Qm = 0.28, Oa = 0.72, h = 
Ho/(100 kms - 1 Mpc -1 ) = 0.7, a primordial power spectrum 
index n =0.96, an amplitude of the mass fluctuations scaled 
to (78 =0.82, and = 0.046. The volume is 32h -1 Mpc on a 
side with 512 3 gas particles and 512 3 dark matter particles. 
This results in a gas particle mass of 4.5 x 10 6 A/q and a dark 
matter particle mass of 2.3 x 10 7 Mq, allowing us to resolve 
galaxies down to stellar masses approaching 10 8 Mq. The 
main simulation output we use for this work is at « = 0.25, 
which is the closest snapshot we have to the typical redshift 
of COS-Halos. 

This simulation includes galactic outflows, which are 
implemented using a Monte Carlo approach and tied to 
the star formation rate (SFR) via M w i n d = 77 x SFR, where 
77 is the mass loading factor. For this work we select two 
wind models to study in detail. The first is our hybrid en¬ 
ergy/momentum driven winds or “ezw” model. This model 
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is our favoured wind model because it most accurately rep¬ 
resents the relevant small-scale physics associated with the 
transition from mom entum-driven winds to energy-drive n 
winds at low masses jHopkinsI 120131 : iMuratov et ail l2015h , 
as well as rep roducing observa tions of the stellar and HI 
mass function (iDave et all2013h . In the ezw model, the wind 
speed v w and mass loading factor r] depend on the galaxy 
velocity dispersion a: 


t) w = 4.29oVfi — 1 + 2.9a 
r / = (j 0 /cr, if a > 75 kms -1 

r / = ((7 0 /(j) 2 , if cr <75 kms -1 


( 1 ) 

( 2 ) 

(3) 


Here, fn = [1.05,2] is the luminosity in units of the Ed¬ 
dington luminosity required to expel gas from a galaxy po¬ 
tential, on = 150 kms -1 , and a is the galaxy’s internal 
velocity dispersion, broadl y constra ined to mat ch IGM en¬ 
richment at high redshift (lOppenheimer fe Davdl2008l l . Par¬ 
ticles are ejected individually, with velocities chosen to lie 
within the range determined by this range of fx,. The range of 
fn valu es is taken fr om ob servations of local starburst galax¬ 
ies by jRupke et al. j2005|) , and is also consistent with obser¬ 
vations by ( Martini 20051 1 that find ft ~ 2. In this model, as 
in its earlier cousin vzw which did not include the steeper 
scaling at low cr, we add an extra kick of 2.9a to simulate 
continuou s pum ping of gas from radiation pressure, as de¬ 
scribed in lOppenheimer fe Pavel |2008t l Q. This simulation 
also includes a prescription for quenching massive galaxies, 
which involves super-heating particles with increasing prob- 
abili ty above a threshold quenching stellar mass as described 
in IDave et all (120131 1 . 

Our second wind model is a constant wind (“cw”) 
model, where the mass loading factors and wind speeds 
are the same for each galaxy, regardless of mass. In the 
cw model, r) = 2 and v w = 680 kms -1 . The latter is 
similar to that used in the Over whelmingly La rge Si mula- 
tions (OWLS) reference model of lSchave et al.l (l2010lj . The 
cw simulat i on us ed here has the same physics desc ribed in 
iFord et ahl (l2013l l , including IWiersma et al] (l2009al l metal¬ 
line cooling, but now with a box size and resolution to match 
our main ezw simulation used here, namely 32/i -1 Mpc and 
2 x 512 3 particles. This allows for consistency with the “ezw” 
model above, and easily resolves galaxies similar in size to 
the smallest ones in the COS-Halos survey. 

The constant wind model does less well reproduc¬ 
ing the ob served measurement s of Civ at high red- 
shift (lOppenheimer fe D ave 2006) , the mass-me tallicity re¬ 
lation ( Finlator et al.l I200S ; IDave et al.l 2011b]}, and the 
evolution of stellar growth in galaxies ( Dave et alj|2011ah . 
Nonetheless; i t via bly reproduces the evolution of global star 
formation ((Springel fe Hernquistll2003l h However, it has not 
yet been tested against the COS-Halos survey at low-z. We 
include this wind model because it serves as a useful contrast 
to the variable mass loading and wind speed wind model 
ezw. Through this comparison we aim to comment on not 
just whether winds are necessary to match CGM observa¬ 
tions, but also what type of winds are necessary. We can also 
gain insight into how baryons cycle in the different wind 


1 Our previous paper listed an incorrect formula for u w ; the one 
shown here is verified to be the one used in our code. 


Table 1. Ion Properties 


Ion 

Ionisation Energy (eV) 

Detection Threshold (A) 

HI 

13.6 

0.2 

Mg 11 

15.04 

0.1 

Sin 

16.35 

0.15 

Cii 

24.38 

0.08 

Sim 

33.49 

0.1 

Si iv 

45.14 

0.1 

Cm 

47.89 

0.25 

Ovi 

138.1 

0.1 


models, how they affect the gas around galaxies, and how 
such processes are manifested in CGM observations. 


2.2 Generating spectra with specexbin 


Once we run our simulations, we use Specexbin to 
calcul ate physical p r operties of the gas along lines of 
sight (lOppenheimer fe Pavel 120061 1. Briefly, Specexbin av¬ 
erages physical properties of the gas along a given sight 
line a nd then u ses look-up tables calculated with CLOUDY 
(|Ferland et al.l ll99sL version 08.00) to find the ionisation 
fraction for the relevant ionic speci e s. We use t he same ver¬ 
sion of_SPECEXBiN_as in [FbrTEjiD (120131 . [2014I ). See Figure 
1 of lFord et al.l (120131 1 for an example of a simulated spec¬ 
trum. This version of Specexbin includes a prescription for 
self-shielding from the ionisation background, which uses a 
density threshold of nH = 0.01cm -3 . Above this density 
threshold we assume that all the hydrogen is fully neutral, 
i.e. it is all Hi, and that all the magnesium is in Mgn. An 
HI column density above ~ 10 17 cm -2 shields against most 
of the 15eV photons that could ionize Mg II to Mg III, but it 
does not shield the 7eV photons that ionize Mgl to Mgn. 

We choose a signal-to-noi se (S/N) value o f 30, for com¬ 
parison with our earlier work (IFord et al .11201 .‘11201 dl l . This 
is higher than the S/N value of the observations, which aver¬ 
ages roughly 12. However, since we will use the same equiv¬ 
alent width detection thresholds as observations, as listed 
in Table [lj our results are not very sensitive to our exact 
choice of S/N. Once the model lines of sight are selected 
(see the next section) and the artificial spectra are gener¬ 
ated using Specexbin, the n we use the Voigt profile fitter 
Aut o VP ( Dave et al.| ll9~9/h to identify the absorption lines 
(see IFord et ~ahl ( 20131 ') for an example) and to fit column 
densities and equivalent widths. 

We choose a large suite of ions to examine for this work 
to get the fullest possible information on the phy sical con ¬ 
ditions of the CGM. As fou nd in earlier work (iWerk et al.l 
120121: IFord et al.ll2013l . 1201414 , a range of ionisation potentials 
are necessary to probe different phases of the CGM. Gen¬ 
erally, low ionisation potential metal lines probe high den¬ 
sity, low temperature regions close to galaxies, while high 
ionisation potential metal lines probe warmer, more diffuse 
gas that can be detected further from galaxies. In addi¬ 
tion to Hll216, we generate spectra for Mgll2796, Sinl260, 
CiU335, Si III1206, Siivl394, Ciii 977, and Ovil032. These 
metal lines are all detected in COS-Halos, and probe a range 
of ionisation energies (see Table [T|). 
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Figure 1. Stellar mass vs. impact parameter for the observed 
COS-Halos survey. Star forming galaxies are separated from pas¬ 
sive galaxies at sSFR = I O ' 11 .V/q fyr. Range bars show the mass 
range of galaxies in the simulations deemed comparable to those 
observed and do not represent errors. 


COS-Halos Data Set 



2.3 Line of sight selection 

We choose our lines of sight (LOS) in our simulations to best 
mimic the LOS in the COS-Halos survey. In COS-Halos, 
galaxies were selected first and the impact parameter to the 
galaxy is the projected distance to the observed quasar. As a 
result, each COS-Halos LOS corresponds to a specific impact 
parameter b around a galaxy of a known stellar mass M„, 
as shown in Figure [l] Moreover, COS-Halos galaxies were 
generally selected to be isolated systems (preferring galaxies 
without photometric coincidences within 1 Mpc, with some 
exceptions), and hence are likely to be central galaxies in 
their halos llTumlinson et al1l2013l . Section 2.5). 

To mimic this process in the simulations, we select 
all the central galaxies in our simulations that are within 
±0.125 dex in stellar mass of the COS-Halos survey galax¬ 
ies, and we generate four LOS at the particular impact pa¬ 
rameter for the selected COS-Halos galaxy, x+b, x-b, y +b, 
and y -b. In this way we produce a suite of model LOS 
for each COS-Halos LOS that matches its unique (6, M*) 
value. For example, for the COS-Halos LOS corresponding 
to M* = 1O 1 o, 2 M0 and b = 18.26 kpc, we select all galaxies 
in our simulation with M* = po 10 ' 075-10 325 , and produce 
four LOS around each at 6 = ±18.26 kpc in x and y. In 
the end, for each of the 44 COS-Halos (6, A/„) LOS, we have 
between 50-200 simulated LOS. The smaller mass galaxies 
tend to have larger numbers of LOS since there are more 
galaxies within ±0.125 dex in M*. 

For consistency with the COS-Halos analysis, we con¬ 
sider all absorption components within ± 600 kms -1 to be 
associated with the galaxy, and sum them for a total equiva¬ 
lent width, unless otherwise stated. W e note tha t this differs 
from our approach in earlier work dFord et al.li2013| ) where 
we only considered components within ± 300 kms -1 to be 
associated with a galaxy. However, as shown in Figures 8 and 
15 of that work, the absorption drops off substantially after 
± 300 kms -1 , so there is little difference between these two 
limits. This same drop-off in absorption after ± 300 kms -1 
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Figure 2. sSFR rate vs. stellar mass for the galaxies in the 
ezw (blue crosses) and cw (red triangles) models that match the 
COS-Halos LOS (black circles). Galaxies with sSFR=0 are plot¬ 
ted as downward facing arrows near the bottom. Over plot ted in 
thick black is the median trend compiled from observations by 
dSpeagle et al.l [201411 . The thin black horizontal line represents 
the dividing line between passive and star-forming galaxies. 


is see n in the COS-Halos data themselves (|Tumlinson et all 
l2013h . 

We note that our simulation naturally produces many 
satellite galaxies around the chosen central galaxy, so in 
some cases our LOS will impact closer to satellites. This will 
presumably contribute to the variance of absorption proper¬ 
ties around a particular galaxy. In this paper, we will focus 
on mean and median absorption statistics without examin¬ 
ing the origin of the scatter, which we leave for future work. 

We also subdivide our sample into passive and star¬ 
forming galaxies (SFGs), using a specific star formation 
rate (sSFR=SFR/M*) criterion of 10 -11 Mg/yr, as done for 
COS-Halos (see Figure [TJ . Figure [2] shows sSFR vs. M, for 
ezw (blue crosses) and cw (red triangles) models, only show¬ 
ing those galaxies selected to mimic the COS-Halos LOS. We 
show as the black line the median trend from a compilation 
of observations of star-forming galaxies by iSpeagle et al.l 
(2014). The observed COS-Halos galaxies are generally ac¬ 
tively star forming for the lower-mass galaxies, but are more 
likely to be passive with increasing stellar mass. Many of the 
massive galaxies in ezw have zero SFR, and those quenched 
galaxies are plotted at the bottom as downward facing ar¬ 
rows. Of the 48 AL* > 10 11 A/q galaxies in the ezw sample, 
all but 9 are quenched, with sSFR=0. Unlike the ezw model, 
cw does not include a quenching prescription, so produces 
no quenched galaxies at high masses. The simulated galaxies 
broadly cover the range of sSFRs seen in COS-Halos, albeit 
that our quenched galaxies have zero SFR. by construction. 

Figure [2] only shows galaxies identified as central galax¬ 
ies within their halos, but at low masses, many of the 
quenched galaxies shown here a re likely former sate llites 
that have been ejected from halos dGabor fc Davel2014l i; our 
quenching prescription is not applied to such small systems, 
but they have been quenched as a satellite likely by gas strip¬ 
ping processes (e.e. IRafieferantsoa et al. I IS). Such satel- 
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lite quenching processes are self-consistently handled in the 
simulations, albeit the version of SPH w e use here do es not 
necessarily ha ndle them optimally ('e.g. lAgertz et ahll2007l : 
IHopkinsI 120131 1. Importantly, we note that the ezw model 
better matches the observed distribution of sSFR vs. M, 
at low masses. This is one reason why we favour our ezw 
model relative to cw, in addition to those arguments dis¬ 
cussed in §[2Tl We also favour the ezw model because it 
more accurately encapsulates the scalings for outflows pre¬ 
dictedfronmecen^Jiiglnresolution individual galaxy simula¬ 
tions iMuratov et al .1120151). and matches many other obser¬ 
vations better than cw I Dave et, al .1120131 ). 


3 COMPARISON TO COS-HALOS 
OBSERVABLES 

We now show a direct comparison between our simulations 
and the COS-Halos survey. This work represents the first 
direct comparison of our simulations against observations of 
absorption lines in targeted LOS at low redshift. Each of the 
following subsections represents a different way of compar¬ 
ing models to observations; all must agree before one can 
claim agreement. We begin by focusing on the strength of 
detected lines, to see if the models match the observed sam¬ 
ple. Then we examine the total absorption along all LOS, as 
well as the incidence rate of all our ions along all LOS. We 
then provide a preliminary glimpse into how ion ratios can 
constrain models, and examine the kinematics of absorbers 
relative to their host galaxies. Matches and mismatches be¬ 
tween all these statistics provide insights into exactly how 
our simulations are succeeding and failing in reproducing 
the COS-Halos data, and they points toward required im¬ 
provements in CGM modelling. 


3.1 Equivalent widths vs. impact parameter 

The most direct absorption line observable is the equivalent 
width. Figure [3] shows the rest equivalent width versus im¬ 
pact parameter for the COS-Halos survey compared to our 
simulations, for our set of probed ions. For this figure we 
focus on detected lines only, and ask if models reproduce the 
same absorption strength as the detected observations. Note 
that the detection thresholds for each ion are given in Ta¬ 
bic |U and shown as solid black horizontal lines to guide the 
eye. We note that for the observations, the detection limits 
quoted are a rough guide, as they can vary vary for indi¬ 
vidual LOS. For the larger sample of LOS in the models, 
however, we use the detection limits quoted in Table [l] for 
all model LOS. Observed points below the detection limits 
are included in this figure for completeness, however in this 
work we only compare observed and model points above the 
detection limits. 

Here we have placed the ions in order of increasing ioni¬ 
sation potential from low to high. The black points are from 
the COS-Halos survey. The circles represent detected values 
while the up and down facing arrows are lower and upper 
limits, respectively. Error bars are included in black for de¬ 
tected values only, errors on limits are not included. In many 
cases the errors are smaller than the size of the symbols. 

The blue crosses and red triangles show the results for 
the ezw and cw models, respectively. These points represent 


the median equivalent width of all LOS with detections above 
the threshold listed in Table [T] at a given impact parameter. 
The bars on the model points span the 16-84% range of EW 
values at that impact parameter, to illustrate the scatter in 
EW at a given impact parameter. To avoid clutter, we only 
show the model range bars for points « every 25 kpc. The 
impact of non-detections will be better illustrated by the 
statistics shown in Figures [4] and [5] 

Both models and observations show a weak trend of de¬ 
creasing equivalent widths with impact parameter. As dis¬ 
cussed in [Ford et alj (l2013l f . the basic driver of this is that 
gas is denser near galaxies. However, the trend is weaker for 
higher ions, and is essentially absent for Ovi, because the 
density gradient is counteracted by the tendency for high 
ions to have increa s ing io nisation fractions in lower-density 
gas. In iFord et all d2013h we found that the ion strength 
dropped with both halo mass (and thus stellar mass) and 
impact parameter. But it is important to note that Figure 
[3] conflates these two trends, owing to the way that we have 
done our LOS selection (0, by matching the observed 
galaxy sample in both stellar mass and impact parameter. 
As shown in Figure [I] LOS at either low or high impact 
parameters can be from either low or high mass galaxies. 
As such, the model points do not decrease as strongly or 
smoot hly with impact parameter, as was seen in lFord et al.l 
d2013lf where we selected narrow bins in halo mass. Addi¬ 
tionally, this degeneracy washes away some trends with ion¬ 
isation potential: the equivalent widths of the lo w ionisation 
species do n ot drop off as steeply as found in For d et ahl 
(|2013Ll20l4l . This will similarly impact Figure [4] The scat¬ 
ter owing to the conflation of these trends is exacerbated for 
the observed points, since the simulation points are averages 
from a larger sample. 

For H I, both simulations are in general agreement with 
the observations across the 0-150 kpc range, showing strong 
absorption in every line of sight. The typical EW is much 
highe r than seen in random lines of sight through the 
IGM dDave et al.ll2010l : lFord et al .1120131 ). showing that there 
is substantial cool gas within these halos. There is little sen¬ 
sitivity to winds, at least as can be discriminated from the 
COS-Halos survey, which is remini scent of the lack of sensi¬ 
tivity to winds in the Lya forest dDave et al.ll201(il '). There 
are, however, a few very strong absorbers with equivalent 
widths above 2A that lie above all the simulation median 
points. These potentially provide an interesting challenge to 
models, hinting at a larger covering fraction of very strong 
Hi than predicted. However, current uncertainties prevent 
any robust constraints emerging from these few systems. 

The metal lines provide a more direct test of pollution 
from galactic outflows. For the low and mid ions, when de¬ 
tections are made, they are generally in the EW range ex¬ 
pected from the models, with a few observed points being 
much higher than the model points (as with HI). Ovi is 
also commonly detected, but the typical equivalent widths 
are generally x2 — 3 higher than predicted; we discuss this 
below. Nonetheless, the broad agreement is a non-trivial suc¬ 
cess of these galactic outflow models in enriching the CGM 
to roughly the observed levels. Though we do not show 
it, a model without winds dramatically fails this basic test 
since it does not yield appreciable metal absorption be yond 
the sm allest impact parameters, as discussed in lFord et al.l 
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Figure 3. Equivalent width (A) vs. impact parameter for observed points (black — circles for measured values, up and down facing 
arrows for lower and upper limits respectively); as well as ezw (blue crosses) and cw (red triangles) models. Model points are the median 
values over LOS with detections above the thresholds given in Table [T] range bars show the 16-84% values of those LOS with detections. 
Ions are ordered according to increasing ionisation energy (right to left, top to bottom). Model points do not smoothly decrease, owing 
to conflation between stellar mass and impact parameter (see Figure [TJ. Errors on observed points are given for detections only (not 
limits) and represent the 2a limits. 


For the low ions (Mgn, Sill, Cn, and Sim), the model 
points tend to lie generally in the middle of the observed 
points. In many cases the scatter in the observed points is 
larger than the 16-84% range bars that the models would 
predict. This could reflect variations in the ionisation condi¬ 
tions or metallicities that are not properly captured within 
our simulations. Recall that we assume a spatially-uniform 
metagalactic ionising background with no local contribution 
from the nearby galaxy and a simple density criterion for 
self-shielding. These assumptions may need to be relaxed if 
we wish to capture detailed variations properly. Doing so 
would require full radiative transfer simulations with much 
higher resolu tion to capture small, dense condensations. As 
discussed in ICrighton et alj (|2015l h if such clouds are in¬ 
deed the dominant contribution to low-ionisation CGM ab¬ 
sorption, they are far from being resolved in any current 
cosmologically-based simulation. 


The mid and high ions tend to show a smaller observed 
scatter than the low ions, more comparable to that seen in 
the simulations. Since low ions tend to a rise closer to galax¬ 
ies and in denser gas liFord et alj|201 j ). this suggests that 
locally-varying density and ionisation conditions, and not 
locally-varying metallicity, is responsible for the larger scat¬ 
ter in the low ions. Alternatively, local variations in metallic¬ 
ity could be much stronger closer to galaxies, which may be 
reasonable if the absorbing gas arises from recent outflows. 

Returning to the discrepancy for Ovi, this is tricky 
to interpret because Ovi ca n arise in bot h p hoto-ionise d 
and collision ally-ionised gas (Qppenheimer fc Pavel 120091 ). 
As shown in lFord et al.l (j20Lj . Figure 6), the median tem¬ 
perature for O VI absorption is ~ 10 4 ' 5 K, which is domi¬ 
nated by photo-ionisation, and indeed the phase space di¬ 
agram in their Figure 5 does not suggest a distinct peak 
of collisionally-ionised O VI absorption at T ~ 10 5 ' 5 K. One 
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way to explain the predicted deficit is thus to invoke a more 
significant hot halo of gas. However, since the typical virial 
temperature of a COS-Halos galaxy is a few times 10 6 K, this 
means that such collisionally-ionised O VI would have to be 
a transient phase through which the gas is cooling. Mod¬ 
elling such non-equilibrium processes is a substantial chal¬ 
lenge for simulations, and may need to be improved if more 
robust and accurate predictions are desired. Alternatively, 
a significantly higher photo-ionising background around 114 
eV than what we have assumed could also result in substan¬ 
tially more Ovi; the metagalactic flux at these energies is 
currently poorly known. 

The details of the wind implementation can strongly 
impact the pha s e of C GM gas. For instance, the models of 
[Hummels et ali l|2013l ) evolved using the adaptive mesh re¬ 
finement hydrodynamic code Enzo show substantially hot¬ 
ter CGM gas, since their models drive winds by super¬ 
heating the ISM. They tended to find that the resulting 
CGM was too hot, and thus also a poor match to Ovi ob¬ 
servations, at least for the most plausible assumptions about 
the energetic input into winds. This highlights that the 
proper treatment of wind propagation, heating, and mixing 
is a challenging problem in current simulations, and shows 
that detailed comparisons to observations can provide inter¬ 
esting constraints on these processes. 

In light of this, it is somewhat surprising that the two 
wind models we consider here are in such good agreement 
with each other, and that these C OS-Halos observat ions can¬ 
not discriminate between them. iFord et~aD (j2013h showed 
large differences in absorption line profiles in every ion be¬ 
tween a cw and a momentum-driven wind model that is 
similar to ezw. Moreover, ezw and cw are well-discriminated 
by other diagnostics generally related to galaxy properties 
liDave et all201 j l. The differences between wind models will 
be slightly more prominent in other CGM statistics that we 
explore later, but for the equivalent widths they are small. 
We will discuss potential reasons for this in §4. 


3.2 Equivalent width per unit redshift 

In the previous section we compared individual absorber 
equivalent widths between the models and the observations 
where detections occurred, and found broad agreement in 
some ions but notable differences in others. That compar¬ 
ison, however, does not characterise the total absorption 
strength; for instance, it may be that absorbers are far more 
or less common in the models compared to the observa¬ 
tions, so that while individual equivalent widths are com¬ 
parable when detected, the overall absorption is discrepant. 
To examine the overall absorption, here we present equiva¬ 
lent widths per unit redshift close to galaxies, including all 
LOS, both detections and non-detections. 

Figure [4] shows the equivalent width per unit redshift, 
dEW/dz, versus impact parameter. We calculate dEW/dz 
by summing the equivalent width (above the detection lim¬ 
its given in Table [TJ in all LOS probed in that impact pa¬ 
rameter bin, and dividing by the total redshift path length 
within ±600 km/s, for all LOS with that impact parameter. 
We bin both observed and model points in 25 kpc bins. The 
observed COS-Halos points are shown in black. For the pur¬ 
poses of this plot, observed points that are upper or lower 
limits are assumed to be at their limit values; the results do 


not change significantly if all upper limits were assumed to 
be zero. The model values are plotted as dotted blue lines for 
the ezw model and as dashed red lines for the cw model. The 
bars on the model points are range bars, showing 16%-84% 
of the model values for all LOS in that impact parameter 
bin (and is likely an overestimate of the actual range at any 
given impact parameter in that bin). The errors on the ob¬ 
served points represent the average errors on all constrained 
detections in the 25 kpc bins. These errors do not include 
errors on LOS that are upper limits, lower limits, or non¬ 
detections; hence the errors shown here are a lower limit on 
the actual error. The errors on the model points have been 
slightly offset horizontally for ease of viewing. 

dEW/dz serves to diagnose whether our simulations are 
producing the observed amount of HI and metal ion absorp¬ 
tion in the CGM around these galaxies. The HI absorption 
is in excellent agreement at all impact parameters. For com¬ 
parison, a random LOS through the universe would typically 
yield dEW/dz«1.24, which is much lower than the CGM 
typical value of several hundreds. Hence the agreement with 
Hi is non-trivial, and shows that there is both the correct 
amount and the correct radial distribution of neutral gas in 
the simulations. 

Turning to the metals, Mg II and C II are well- 
reproduced in the model, with the observed points being 
within the 16 — 84% range of the simulation predictions. 
The models produce drastically too little Si II and Si III; the 
models are off by as much as an order or magnitude. Con¬ 
versely, the models produce too much Siiv at low impact 
parameters, even outside the model ranges. The simplest 
interpretation is that the models do not have as much cold, 
dense gas in the CGM as observed, while having an over¬ 
abundance of warmer gas. A physical interpretation might 
be that our assumption of optically-thin gas fails to capture 
the dense gas phase structure in the CGM. Qualitatively, 
self-shielding would be expected to increase Si n and reduce 
Si IV. However, we note that substantial self-shielding would 
greatly increase the amount of H I and C II, which at present 
are in good agreement with the observations. Clearly there is 
some disagreement in the phase structure predicted in the 
models compared to that observed, but the nature of this 
remains difficult to ascertain. Small features in the CGM 
which may give rise to absorption, such as interface regions 
and small, possible multi-phase substructures, may be below 
the resolution limit of the models presented here. 

For the mid-ion C ill, the models produce slightly too 
little absorption but generally the observed points are within 
the model range values. For Ovi, however, the disagreement 
is much more stark. While the difference in the strength of 
detected lines (as shown in Figure El) between the models 
and the observations is roughly a factor of 2, the difference 
in the total absorption over all LOS is more like a factor 
of 5. This difficu l ty in p roducing enough Ovi to match the 
iTumlinson et al] (1201 If) observations was a lso found in the 
Enzo simulations of Hummels et al.l (l2013fl . 

The agreement between the models and the observa¬ 
tions for H I, Mg II, C II, and to a lesser extent C ill is gener¬ 
ally quite good given that the models have not been tuned 
in any way to reproduce these data, and suggests that the 
amount of CGM enrichment in both these simulations is 
roughly in accord with observations. However, the mismatch 
in the silicon ions, as well as Ovi, is quite large. This is a 
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Figure 4. Summed equivalent width divided by total path length (dEW/dz) over all LOS at a given impact parameter, as a function of 
impact parameter, binned in 25 kpc bins. Black points are from the COS-Halos survey, with error bars averaged from well-constrained 
detections only, no limits or non-detections. Blue dotted lines are for the ezw model and red dashed lines for the cw model. Both models 
have range bars showing 16%-84% of all model LOS in a given impact parameter bin with detections above the thresholds given in Table 
|T| The errors on the model points have been horizontally offset slightly for easier viewing. 


significant failure, given both the degree of the mismatch, 
and the high quality of the observed Si ill and Ovi data. 
The similarity in the predictions of the two wind models is 
striking, particularly since the amount of CGM metals and 
the phase space distribu tion of the gas in these two models 
is rather different dFord et alJhoidlf . as we will discuss later. 

The trend in HI is quite flat with impact parameter. 
As we discussed in the previous section, this partly owes to 
the fact that a given impact parameter bin has a substantial 
range of stellar masses. Some of this likely owes to satura¬ 
tion, as HI at these equivalent widths are very much on the 
linear part of the curve of growth and hence fairly insensitive 
to variations in column density. Nonetheless, the trend in HI 
with b is even much flatter than in any of the low-ionisation 
metals. At face value, this suggests that the dropoff in metal 
ions at large radii owes either to gas becoming more ionised 
at higher impact parameter or to a dropoff in metallicity, 
rather than a dropoff in the amount of neutral gas. 

In summary, the predicted total equivalent width per 
unit redshift shows broad agreement with the observations, 
but with various interesting discrepancies. Hi, Mgn, Cn, 
and C ill are in fair agreement with the observations, and 
Si II, Si ill, and O VI are significantly underproduced. It is dif¬ 
ficult to identify a single physical variation that could even 
qualitatively reconcile all these discrepancies while preserv¬ 


ing the good agreements, suggesting that further work is 
needed to discern their detailed cause. Radial dependencies 
are stronger for lower ions, except for HI that shows almost 
no decline with radius possibly owing to saturation. 


3.3 Covering Fractions 


The covering fraction quantifies how patchy the CGM ab¬ 
sorption is in a particular ion. Observations and simulations 
suggest that the CGM has a complex, multi-phase struc¬ 
ture, potentially with s mall clouds embedded in warmer 
gas llTri pp et al. 20081: Kacprzak et alJ l2012j ; iFord et ah! 
120141 ; Crighton et al] 2015ll . The size and fraction of ma- 
terial in these variou s phases remains poorly constrained. 

I Churchill et al.l d2003if . based on observed Mgn absorption 
around galaxies, finds that this absorption generally arises 
in small clouds, i.e. 2 ^ uh ^ 20 cm -3 , line-of-sight sizes 
1 ^ D ^ 25 pc, and m asses between 10 and 1000 M/Mq. In 
contrast. IWerk et al.l J2014I ). from low-ionisation absorbers 
in COS-Halos, find a median cloud mass of 10 7,6 Mq. Cov¬ 
ering fractions thus provide another test of models, and 
potentially a very stringent one if the cloud sizes are well 
below what can be resolved in our (or most) simulations 
dCrighton et alj[2015h . 

Figure [5] shows the covering fraction as a function of 
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Figure 5. Covering fraction of COS-Halos survey data (black solid lines are values, black dotted lines are errors), compared to values 
from the ezw (blue crosses) and cw (red triangles) models. Observed values are in 50 kpc bins, and model points are not to show the 
downward trend with impact parameter. 


impact parameter for our set of ions. As in Figure]]!] we plot 
ions from low to high ionisation potential. We define a cov¬ 
ering fraction as the fraction of LOS with a detection above 
a given ion-specific EW threshold, corresponding approxi¬ 
mately to the COS-Halos survey detection limits, as listed 
in Table]]] The observed points are binned from 0-50 kpc,50- 
100 kpc, and 100-160 kpc. The solid black line represents the 
median covering fraction, and the dotted black lines repre¬ 
sent the uncertainties (Wilson intervals). The models are 
plotted as individual points (crosses for ezw, triangles for 
cw, as before), because in the simulations we have multi¬ 
ple LOS at each impact parameter and are able to define a 
unique covering fraction at each impact parameter. In this 
manner one can get a sense of the variability in covering 
fraction as a function of impact parameter, albeit with the 
same degeneracy in stellar mass noted earlier. 

For Hi, the simulations almost always have unity cov¬ 
ering fraction above its detection threshold of 0.2A. The 
covering fraction in the COS-Halos survey is also very close 
to unity, but formally slightly bel ow. This high HI frac¬ 
tion was also emphasised in lProchaska et all (1201 11. and is 
independent of both galaxy stellar mass and galactic star 
formation rate (not shown here), for the range of these val¬ 
ues covered by the COS-Halos survey. In fact for all ions 
shown here, the covering fraction depends more strongly on 
the impact parameter than either the star formation rate or 
the stellar mass. The high covering fraction of Hi is noted 


in the COS-H alos survey for both star-forming and passive 
galaxies dThom et al.,' .20121 ). We note that we are able to 
reproduce these high HI covering fractions in passive sys¬ 
tems with our heuristic quenching prescription (present in 
the ezw model but not the cw model). 

Further exploring Figure [5] we see that for low ions 
(Mg II, Si II, C II, and Si ill) the models generally are either 
at, or slightly above, the observations at impact parameters 
of 0-50 kpc but drop below the observations at larger im¬ 
pact parameters. For the mid ion Si IV, the model values are 
below the observed values in the 0-50 kpc bin, but within 
the errors. At higher impact parameters the model covering 
fractions drop below the observed errors. For C ill, there are 
model points below the observed errors at all impact param¬ 
eters. For the high ion O VI however, the model covering frac¬ 
tions are already too low at all impact parameters, with the 
match worsening at higher impact parameters. We note that 
the covering fraction of the high ion O VI does not decline 
as steeply with impact parameter as the lower ionisation 
species, either in the observations or in the models. We also 
emphasise the similarity of the ezw and cw model points. 
While there is scatter in the model points, there is not a 
single impact parameter bin in a single ion where one model 
is consistently different from the other. Rather, both models 
match (or mismatch) the observations in similar ways. 

The mismatch between the models and the observations 
has three possible origins: there could be not enough metals 
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in the simulation, the metals could be in other ionisation 
states not shown here, or the metals could be in the wrong 
spatial distribution, which could happen if the metals were 
not mixed properly. We note that in both Figure [4] and [5] the 
radial gradients of Sill, Sim, Cm, and Ovi are steeper in 
the models than in the observations. This suggests the mod¬ 
els do not produce enough metals in the outskirts of halos 
to match the COS-Halo observations. It is also possible that 
some of the metals are not in the correct ionisation states: 
note our overproduction of Si IV at small impact parameter 
but an underprediction of Sill and Sim. We also slightly 
overpredict C II but underpredict C III at low impact param¬ 
eters, although the discrepancies are generally within the 
error bars of the observations. It is likely a combination of 
these factors that account for the mismatch, as we discuss 
in further detail in i)3.5l and i )4.2l 

In our simulations, we do not mix metals into CGM gas 
from the particles that are ejected from the ISM, and hence 
it is perhaps not surprising that there is too little metal mix¬ 
ing in our models. It is not obvious how much mixing should 
occur; again, it is challenging to model mixing of outflowing 
gas properly. Nonetheless, it is interesting that the covering 
fraction can potentially provide a good constraint on this 
crucial physical process affecting how the outflows and the 
CGM interact. 

3.4 Dependence on the Specific Star Formation 
Rate 

A key feature of COS-Halos is that it includes a significant 
sample of passive galaxies, and hence can examine trends 
of CGM absorption with star formation rate. Current mod¬ 
els broadly suggest that passive g alaxie s are s urrounded by 
hot gaseous halos (e.g. lOabor & Dava 120141 1. unlike star- 
forming galaxies t hat are located in “cold mode” halos (e.g. 
iKeres et alJ l2005f l. Hence one might expect trends in the 
strength of various ions as a function of specific star forma¬ 
tion rate. 

We examine the dependencies on sSFR by investigat¬ 
ing equivalent widths versu s impact parame t er, sep arated 
by sSFR. As noted earlier, iTumlinson et al] (120111 ) found 
that blue galaxies generally have higher equivalent widths 
of Ovi than red galaxies, r egardless of impact parameter 
or galaxy mass. In contrast, iThom et al.l J2012ll found that 
the HI properties in the CGM of red and b lue ga laxies are 
only marginally different, and IWerk et al.l ( 201 similarly 
found that low-ionisation absorption is also not greatly dis¬ 
similar. We use our ezw model, which includes our heuristic 
quenching prescription, to examine how well our simulation 
reproduces these observations. 

In Figure [6] we divide our model LOS - which are ini¬ 
tially selected based only on impact parameter and stellar 
mass - into LOS around high- and low-sSFR galaxies (above 
and below 10 _11 MQ/yr). To focus the discussion, we exam¬ 
ine only HI, our lowest metal ion Mgn, and our highest ion 
O VI; the trends for intermediate ions are interpolatable from 
these. 

We are able to reproduce t he observed high HI covering 
fractions in passive systems dThom et al.l 12013 1 despite our 
heuristic quenching prescription. In fact there is little dif¬ 
ference in the equivalent width of H1 in model red and blue 
galaxies (top panel), which is consistent with what is seen 


in the observations (blue stars and red crosses). Hence our 
quenching model does not remove HI from halo gas, even 
though it adds substantial energy to the halo by heating 
the ISM gas to temperatures that are well in excess of the 
virial temperature. We are not claiming that this lends phys¬ 
ical validity to our simple quenching model, but rather that 
such observations could provide interesting constra ints on 
more physically-motivated models of quenching (e.g. ICaboi] 
120131) . A similar trend is seen in Mgn (middle panel), where 
our simulated red and blue galaxies have similar absorption 
strengths. In the observations, if anything there seems to 
be more Mg 11 around red galaxies. Our results suggest that 
whatever quenches the central galaxy in more massive halos 
should have a minimal effect on the cool CGM gas that gives 
rise to HI and Mgn. 

For Ovi, the situation is different. The observed 
strength around bl ue galaxies a re higher than around red 
galaxies (jTumlinson et all 1201 il l, three of which are upper 
limits off the bottom of the plot. In contrast, this difference 
is not present in the simulations: there is no clear difference 
between the equivalent widths of O VI in red and blue galax¬ 
ies. Hence, the underprediction of Ovi in the models seen 
earlier is specific to blue galaxies, while the models greatly 
overpredict Ovi for passive galaxies. 

This shows that, despite ejecting gas at high tempera¬ 
tures from the ISM, our quenching model probably does not 
properly heat the surrounding CGM gas to prevent cool¬ 
ing through the O VI temperature regime. Our implemented 
quenching model does not prevent accretion by keeping sur- 
roundi ng gas hot , but o nly ej ects the gas once it gets in the 
galaxy. iGabor fc Pavel (l2014l j argued that keeping the hot 
halo gas near the virial temperature, well above where it 
would be in the Ovi phase, provides a successful model to 
produce the quenched galaxy population as observed. Im¬ 
plementing such a model into our simulations would likely 
have quite a different impact on O VI absorption, as it would 
explicitly prevent gas from cooling. It remains to be seen 
whether such a model would concurrently over-suppress H1 
and low ion metal absorption; we leave this investigation for 
future work. For now, we conclude that our implemented 
form of quenching is inconsistent with CGM metal obser¬ 
vations. Clearly, matching both the low ions and O VI in 
the CGM of passive galaxies represents an important con¬ 
straint on the thermodynamics of quenching. Although pas¬ 
sive galaxies in our simulations are not currently ejecting 
winds, we show in §4.3 that the CGM absorption in our 
models is dominated by gas that was ejected well before the 
epoch of observation, largely decoupling the CGM predic¬ 
tions from the current star formation rate. 

3.5 Ion Ratios 

Ion ratios provide a complementary test for the physical 
conditions of the absorbing gas. Ratios between metal ions 
and hydrogen provide an estimate of the metallicity in the 
gas traced by that ion, while ratios between ions of the same 
element provide a constraint on the ionisation conditions of 
the gas. Given our current analysis, we can only study ratios 
of ion equivalent widths summed over 600 kni/s intervals, 
which likely erases some of the detailed information about 
physical conditions through the CGM; trying to subdivide 
these into components introduces great sensitivity to the 
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Figure 6. Equivalent width vs. impact parameter of H I, Mg II, and O VI in the ezw model. Red model galaxies are shown as red diamonds, 
blue model galaxies as blue squares, red observed galaxies as red x’s, and blue observed galaxies as blue asterisks. Red and blue are 
separated by a sSFR = 10 —11 Mq/ yr. Range bars show the 16-84% values of the model sightlines at fixed impact parameter, as in Figure 
[3j but are larger here because the sample has been divided. 


exact noise level and fitting procedure, which is beyond the 
scope of this work to consider. Nonetheless, we will see that 
summed ion ratios still provide interesting constraints on 
CGM properties in our simulations. 

Selected ion ratios are shown in Figure [3 As before, 
COS-Halos survey data is plotted in black, and model points 
(ezw: blue crosses, cw: red triangles) represent the median 
value of the ratio at each impact parameter. For the observed 
points, downward facing arrows indicate that the ratio is an 
upper limit (meaning an upper limit divided by an actual 
value, or an actual value divided by a lower limit), upward 
facing arrows indicate the ratio is a lower limit, and circles 
indicate an actual value. We do not plot values that are 
poorly constrained, such as a limit divided by a limit. For 
the model points, coloured downward arrows show where 
the ratio is zero, as little or none of the top ion is present. 
In the left panels, we plot the log of the ratio of equivalent 
widths of a low and high ion relative to Hi, as a probe of 
the metallicity of the gas traced by that range of ionisation 
potentials. In the right panels, we show the ratio of two Si 
lines, as well as a ratio of a low (Mgn) to high (Ovi) line. 

The low ion ratio (Mgn/Hl) shows generally good 
agreement with the observations, albeit with a scatter in 
the observations that is much larger than in the models. 
The scatter might arise owing to local variations in ionisa¬ 
tion conditions. There is a slight hint that the models pro¬ 


duce slightly too high a ratio, particularly at large impact 
parameter, indicating that the metallicity in the low ionisa¬ 
tion gas might be a bit too high. There is also a non-trivial 
offset between the wind models, ezw is slightly higher than 
cw, although the scatter in the observations is far too large 
to discriminate between the models. Still, this is one of the 
few aspects in which any statistically significant difference 
can be seen between the ezw and the cw models. 

The predicted Ovi/Hl is considerably lower than most 
of the observed points. The discrepancy is comparable to 
that seen in the individual EW measurements of O VI vs. im¬ 
pact parameter in Figure [3] These discrepancies suggest that 
the hot diffuse phase in our models is either insufficiently en¬ 
riched with metals or has an incorrect temperature structure 
that shifts oxygen to higher or lower ionisation states. These 
trends hold true for both wind models (though ezw produces 
slightly higher ratios than cw), indicating that it is not spe¬ 
cific to which wind model we use, but rather may be indica¬ 
tive of an overall failure in the way we implement winds, 
namely by ejecting cool, unmixed ISM particles. Ejecting 
hotter or over-enriched gas that mixes would qualitatively 
be expected to move towards better agreement with the ob¬ 
servations. Nonetheless, we emphasise that the results are 
within a factor of « 3 of the observed ratios, which is already 
a non-trivial success. These ion ratio comparisons further re¬ 
fine and highlight issues with our models, and demonstrate 
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Figure 7. Comparison between the log of the ratio of equivalent width in the COS-Halos survey (black symbols: circles for constrained 
values, downward arrows for upper limits, upward arrows for lower limits), as well as for the ezw model (blue crosses) and the cw model 
(red triangles). Coloured downward triangles show where the model ratios are zero. Left panels show various metals vs. Hi, probing 
metallicity. Right panels show metal ratios, probing physical conditions. 


the ability of the COS-Halos survey to provide important 
insights into galactic outflow and ionisation processes in the 
CGM. 

The right panels of Figure 0 plot different metal ions 
against each other, mitigating the dependence on metallicity 
and allowing us to test physical conditions of the gas more 
directly. However, the observational data are more limited, 
as it is not common to find absorbers with multiple ions 
detected. The ratio of Siii/Sim at low impact parameters 
shows an overall trend of increasing ratio with impact pa¬ 
rameter, and both models and are in good agreement with 
the observations from ~ 0-50 kpc. Beyond 100 kpc, the ezw 
model has a higher ratio than the observations, indicating 
cooler gas than observed. Meanwhile, the cw model shows 
a very low Si Il/Si ill ratios, typically having essentially no 
Sill and hence dropping off this plot at these large impact 
parameters (shown as coloured downward arrows). The scat¬ 
ter is less than in the ratios versus Hi, suggesting that the 
ionisation conditions giving rise to particular ions are more 
uniform in the CGM. 

The Mgn/Ovi ratio shows a large scatter in the ob¬ 
servations, indicating that these two ions arise in different 
phases that do not trace each other well. There is, perhaps 
surprisingly, no clear trend with impact parameter, as one 
might expect for a simple model in which the low ion arises 
only near the central galaxy and the high ion arises through¬ 


out the halo dFord et al.ll2013h . This suggests that the Mgn 
is occurring throughout the halo, perhaps in satellite galax¬ 
ies or cold dense clouds interspersed with warmer diffuse 
gas. Meanwhile, the models show a rising trend at small 
impact parameter that flattens at larger impact parameter. 
The median ratio does not show a large scatter, suggesting 
that Mg II and Ovi are appearing together more often in 
these models than in observations. 

In summary, by examining ion ratios we can isolate 
trends with metallicity and ionisation state in the CGM as 
traced by various ions. While there is broad agreement be¬ 
tween the models and the observations for HI and some low 
ionisation metal species, there is also significant disagree¬ 
ment for higher potential ions, particularly Si ill and O VI. 
Clear trends emerge that suggest that the metallicity within 
certain phases may be too high (for warm diffuse gas traced 
by mid ions) or too low (for hotter diffuse gas traced by high 
ions). Ion ratios generally show moderate agreement with 
the observations, albeit with smaller scatter and perhaps 
some discrepancies at large impact parameter (> 100 kpc) 
where the low ions disappear in the observations. The dif¬ 
ferences between our two wind models are small compared 
to the differences between the models and the observations, 
which may owe to the similarly kinetic and cool ejection of 
material. Nonetheless, the differences between wind models 
seen here are not seen in most of the other statistics we 










Cycling in Low-Z CGM 13 


have considered so far. Ion ratios thus provide new and in¬ 
teresting challenges to models of CGM gas, and obtaining a 
larger and higher-S/N sample could potentially provide the 
best constraints on the physical conditions in the CGM. 


3.6 Kinematics 

With spectroscopic redshifts for all its galaxies, the COS- 
Halos survey provides complementary information on the 
kinematics of absorbing gas relative to the systemic velocity 
of the host galaxy. In principle, kinematics could distinguish 
outflowing gas fr om inflowing, but in practice this division 
is not so clean (IFord et al.ll2014f ). Nonetheless, the kinemat¬ 
ics provide an additional constraint on CGM models. In this 
analysis we will no longer sum absorption within ±600 km/s 
of the galaxy, but instead consider individual components to 
elucidate the kinematics. We note that this introduces ad¬ 
ditional sensitivity to how our profile fitter decomposes the 
absorption into components, but here we really focus only on 
kinematics and not line strengths, hence such a sensitivity 
should not greatly impact our results. 

In Figure [5] we quantify from where in velocity space 
the absorption originates in the observed data, and how 
this compares to the model predictions. We focus on the 
zw model; the cw model gives very similar results. To do 
so, we lay out a coloured “chequerboard”, where we bin 
our model data into 25 kpc spatial bins and 100 kms -1 
velocity bins. For each of those bins, we calculate the to¬ 
tal column density for all LOS with impact parameters in 
that bin, for all absorption within ± 600 kms -1 . We shade 
each block by the fraction of the column density in that 
impact parameter range that comes from a given velocity 
interval. Dark blue blocks indicate that more than 30% of 
the total column between ± 600 kms -1 comes from that 
given velocity interval at that impact parameter, light blue 
blocks indicate 20-30%, green 10-20%, and orange that less 
than 10% comes from that velocity interval. We note that 
these percentages are for that specific impact parameter bin 
only, and not for the entire 0-150 kpc range, so each vertical 
column of shaded blocks is independent. We also plot the 
location of each individual observed absorber (not summed 
over velocity, as elsewhere in this work) as open circles whose 
sizes indicate the strength of the absorption. For the metals, 
the large circles indicate log(lV) > 14.5, the medium circles 
log (N) = 14.0 — 14.5, and the small circles log (N) < 14.0. 
For Hi, the scale is shifted up by one in the log. 

Generally, all the observed absorbers lie in the range 
where the models predict a significant fraction of the ab¬ 
sorption should lie. The observed absorption also clusters 
strongly around the systemic velocities of the host galaxies, 
which is in line with the model predictions. The exception 
to this is that at very low impact parameters there is an 
excess of absorbers with positive velocities, and for HI and 
C ill that absorption is stronger than what is expected from 
the models. This owes to a single sightline that happens to 
have strong absorption in multiple ions and lies at a positive 
velocity with respect to its host galaxy. 

Overall, both the simulations and the observations sug¬ 
gest that most of the absorption occurs roughly within the 
virial velocity of the galaxy’s halo. Only absorbers with the 
smallest impact parameters show any sign of an outflow, 
and it is difficult to assess its significance because it only 


occurs in one system (with several absorbers within it), and 
there could be systematic uncertainties identifying the kine¬ 
matic centre of any given system. Otherwise, the simulations 
broadly match the observed kinematics, and do not indi¬ 
cate an abundance of strongly outflowing gas at the present 
epoch, in g eneral agreement wi t h the outside-in enrichment 
scenario o f lOppenheimer et alj (l2012lb 


4 COMPARING THE CGM IN WIND MODELS 

In general, we have found rather modest differences in obser¬ 
vational tracers of absorption as probed by COS-Halos be¬ 
tween our two wind models. Both models give similar predic¬ 
tions for the equivalent width of detected lines, total equiv¬ 
alent width, and covering fraction of HI and some of the 
low ionisation metal lines, which broadly match data. Both 
models fail, most notably, to reproduce the total equivalent 
widths of Sim and O vi, and do so in similar ways. While we 
favour the ezw model for independent reasons, in the CGM 
observations presented here it typically does no better or 
worse than the constant wind model in matching the obser- 
vations. This is somewhat curious, since earlier w ork (e.g. 
lOppenheimer fe DaveH2008l : IFord et al. 1120131 . 12014) showed 
significant observable differences between these wind mod¬ 
els in the physical conditions and enrichment of CGM gas. 
We examine the reasons and implications of this surprising 
agreement here. 


4.1 Halo Mass-Stellar Mass Relationship 


A key difference relative to our previous works is that here 
we are obtaining a galaxy sample that is matched in stellar 
mass, not halo mass. In our previous work we chose to match 
ezw and cw in halo mass since we reasoned that the CGM 
is more tied to the halo, e.g. the chance of metals escaping 
the CGM is tied to the full halo potential not just that of 
the galaxy, and the total amount of CGM baryons should 
scale more with the halo than with the galaxy. However, 
COS-Halos has chosen galaxies based on their stellar masses, 
and hence in this observational comparison we have focused 
on matching this quantity. As we argue here, it turns out 
that the “better than expected” agreement between the two 
models presented here owes primarily to this choice, which 
has interesting implications. 

Figure [9] shows the stellar mass-halo mass (SMHM) 
relation for our two wind models, with a running median 
for each (blue dotted for ezw, red dashed for cw). We in¬ 
clude galaxies down to our stellar mass resolution limit of 
1.4 x 10 8 Mq (32 gas particle masses). Solid black lines de¬ 
lineate the stellar mass range of COS-Halos. The broken 
coloured lines show results from abu ndance matchin g anal¬ 
yses for the closest available red shift: BehroozieMilJ (|2013l l 
as the green dash-dot line, and lMoster et afl ~( 20131 ') as the 
purple dash-dot-dot line. The ezw model is clearly a bet¬ 
ter match. This simply reflects the fact that ezw (or its 
cousin, vzw that does not have the energy-driven scalings at 
low masses or quenching prescription at high masses) pro- 
vid es a better match to the ga l axy stellar mass fu nction than 
cw (lOppenheimer et ai1l20ld : iDave et alJl2011ah . 

The stellar mass range relevant to this wor k is 
1q9.5— h. 5 ( S ee Figure [9]), which from iBehroozi et alJ 
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Figure 8. Distribution of column density along the line of sight. Ezw model data has been binned both in impact parameter and in 
velocity space (cw, not shown here, is similar). Dark purple blocks show that 30% or more of the column density between -600 and +600 
is concentrated in that velocity bin. Blue shows 20-30%, green between 10-20%, and orange less than 10%. Overplotted as black circles 
are the COS-Halos observed points, scaled by the absorber strength. For the metals, large circles show log(TV) > 14.5, medium circles 
show \og(N) = 14.0 — 14.5, and small circles show log (TV) < 14.0. For Hi, the scale is increased by 1: large circles show \og(N) > 15.5, 
medium circles show log(Af) = 15.0 — 15.5, and small circles show log(Af) < 15.0. 


J 2013 I 1 corresponds to a halo mass range of 10 11 ' 4 15+ Mg. 
In much of this range, the cw model’s SMHM relation is 
lower than ezw’s. This mainly occurs because the cw model 
ejects material at higher velocities in this mass range, which 
reduces the amount of wind recycling that dom inates the 
inflow at these masses (jOppenheimer et al.ll201Ci l. 

I 11 the present work, we compare the wind models at 
fixed stellar mass. Hence for a given M„, the cw model will 
predict a higher Mhaio- We can use this to qualitatively rec¬ 
oncile whv lFord et al.l (120131 1 found large differen ces between 
the tw o wind models. For example, Figure 9 of lFord et al.1 
d2013h shows that for a momentum-driven wind model (sim¬ 
ilar to the ezw model used here), higher mass halos have 
higher column density and hence EWs at a fixed impact pa¬ 
rameter. Figure 14 of that same work shows that for a fixed 
halo mass, the constant wind model has a fewer column den¬ 
sity at fixed impact parameter. Therefore, at fixed stellar 
mass, given that the constant wind model has higher halo 
masses than ezw, this would result in a higher column den¬ 
sity at fixed impact parameter. The differences in absorption 
strength owing to the wind model are thus offset by varia¬ 
tions owing to halo mass. As a result, for the observables 
presented above, the two wind models match much better 
at fixed stellar mass than at fixed halo mass. This then prop¬ 


agates into all the various other statistics, such as covering 
fractions and dEW/dz. 


The interesting physical implication of this is that our 
models predict the CGM metal absorption properties are 
more tied to the stellar mass than to the halo mass. On 
the one hand, this is surprising since the CGM represents 
halo gas, which one might think would be more closely tied 
to the halo mass. That being said, the metals in the CGM 
are formed by the stars in the galaxies and then ejected, 
so the closer r elatio nship is perhaps not une xpected. In 
iFinlator et al.l (120081 4 and iDave et all (l2011bh it was ar¬ 
gued that the metals retained within the ISM, as traced 
by the mass-metallicity relation, are strongly governed by 
outflows, whose properties must be tied most closely to 
stellar mass in order to give rise to a t ight s tellar mass- 
metall i city relation (e.g. Finla tor et al] 120081 : IDave et al.l 
l2011bl : IPeeples fe Shankai 201 ill . Evidently, in these sim¬ 
ulations, the metals deposited into the CGM are likewise 
strongly governed by outflows that are more closely tied to 
the stellar mass. 


We emphasise that this similarity in predictions for ezw 
and cw does not imply that current CGM observations are 
unable to rule out any outflow models. While we have varied 
the free parameters in our particular implementation of out- 
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log (M hal0 ) 

Figure 9. Stellar mass-halo mass relationship for central galaxies in the ezw (blue dots, averaged as the dotted line) and cw (red dots, 
averaged as the dashed line) at z=0.25. The stellar m ass range of th e COS -Halos data set is bounded by solid black lines. R esults from 
abundance matching are shown: green dash-dot line llBehroozi et a.l.ll20 1 -'ll and purple dash-dot-dot line llMoster et al.l[20l1 ). 


flows, there are many other implementations of outflows in 
simulations that may result in quite different CGM proper¬ 
ties. These include model s that over pressurise the ISM by 
local supernova heating (IHummels et alJ l201.il : IShen et al.1 
120131 : iRoskar et al.l I20l4) , drive winds thro ugh a combina - 
tion of radiation pressure and supernovae (lHopkina 120131) . 
and include early _stellar fee dback driven by stellar winds 
and radiation (IStinson et al.ll2013l ). Compared to these, our 
simulations make some rather extreme assumptions, namely 
that the ejected particles are not initially heated (though 
they can be shock heated by interactions with the ambient 
halo gas), and that our metals do not mix with the ambient 
gas. It appears that these assumptions may be more critical 
to setting the CGM enrichment and phase structure than the 
details of the assum ed wind speed and mass loading factor. 
For instance, iHummels et al.l (|2013h found that their AMR- 
based models did not fare well against CGM observations, 
except with rather extreme assumptions about outflow en¬ 
ergetics. We look forward to detailed comparisons of other 
simulations, similar to what we have done here, in order to 
better assess how the physics driving outflows manifests in 
CGM absorption line data. 


4.2 Baryonic Fractions within Halos 

The two wind models considered here predict similar CGM 
absorption line properties, despite differences in the amount 
and physical state of CGM gas. We examine those differences 
here. 

Figure [lO] shows the cumulative enclosed baryonic frac¬ 
tion (top panel) and metal fraction (lower panel) within 
halos as a function of radius, for our COS-Halos matched 
sample of galaxies. The upper panel shows the total bary¬ 


onic mass fraction divided by the cosmic baryon fraction (i.e. 
B(,/O m ). Unity indicates a cosmic mean baryon fraction, and 
shows approximately what models would predict in the ab¬ 
sence of outflows (in detail, it wou ld be slightly higher than 
unity at these masses; see e.g. (Day il l2009l ). The solid black 
line shows the total fraction of baryons (stars+gas) in ha¬ 
los for the ezw model, as a function of enclosed radius from 
the galaxy. Even at the virial radius, this is less than unity, 
showing the efficiency of outflows in driving out baryons. 
The dotted black line shows the baryon fraction for the cw 
model, which is substantially lower than that of the ezw 
model. Owing to higher wind speeds in this halo mass range, 
the constant wind model carries more of its baryons out of 
halos into the IGM than ezw. Both models push baryons 
out preferentially from the central regions, where most of 
the star formation (and hence outflow generation) occurs. 
This leaves a large component of gas-phase baryons in the 
CGM. 

For the ezw model, CGM gas (hot+cool phases) makes 
up 65% of all baryons inside the halo. This is broadly con¬ 
sistent with estimates of ba ryonic mass derived from ob¬ 
servations (|Werk et al.ll2014l ). The cw model predicts many 
fewer baryons, particularly in the cool phase. However, it is 
important to note that in converting observable quantities 
into estimates of mass, assumptions must be made about the 
physical conditions in order to obtain ionisation corrections. 
The consistency of ezw here, then, may simply reflect the 
similarity in phy sical con ditions in the ezw model to those 
assumed in (iWerk et al.ll2014l ). An assumption of a higher 
proportion of hot gas, as present in the cw model, would 
alter the conversion of observable quantities to mass. Hence 
this should not be regarded as a success of ezw relative to 
cw, but may provide guidance on what assumptions should 
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Figure 10. Baryonic content (top panel) and metal content (lower panel) of our COS-Halos matched galaxies for ezw (solid) and cw 
(dotted) wind models as a function of enclosed radius. The black lines show the total baryon fraction for the ezw and cw models. Dark 
blue, red, light blue, and orange lines show the fractions of ISM (star-forming) gas, hot gas (> 10 5 K), cool gas (< 10 5 K), and stars, 
respectively. 


be made in order to obtain CGM mass estimates from ab¬ 
sorption line data. 

The ISM of galaxies in the cw model contains an order 
of magnitude less baryons than in the ezw model, and about 
one-third as many stars. The reasons for this are two-fold: 
more gas is driven out of the ISM, making less available for 
star formation and strong outflows can actually heat the sur¬ 
rounding region thereby preventing the accretion o f fresh gas 
that would otherwise replenish the loss llOppenheimer et ali 


l20ld : lvan de Voort et al • 2on|). The radial profile of the ISM 
and stellar fractions are both quite flat, since these phases 
are concentrated in the galaxy’s central regions. 

More directly relevant to CGM absorption lines is the 
cool CGM baryon fraction, which is about three times higher 
in the ezw model than in the cw model and shows a modest 
increase with radius. In contrast, the hot (T > 10 S K) baryon 
fraction is not very different in the two models, and increases 
strongly with radius. The fact that the HI equivalent widths 
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are much less than a factor of three different in these two 
models (see Figure EB is, therefore, curious and suggests 
that the global CGM cool gas content is not straightfor¬ 
wardly traced by HI absorption. However, this again may 
owe largely to saturation effects in HI. 

The lower panel of Figure flOl shows the metal fraction 
in various halo gas phases. Here we have plotted the metal 
mass fraction enclosed within a radius R/R v i r , divided by 
the cosmic baryon fraction (i.e. flb/f! m ) for ease of compar¬ 
ison with the upper panel. The constant wind model has a 
lower enclosed metal fraction, since more of the metals have 
been driven out of the halo by outflows and since it made 
less stars. Of the metals that are present in the cw model, 
more are hot. For a sense of the metal content derived from 
observations o f CO S-Halos galaxies, we refer the reader to 
IPeeples et al . (2014), but again we note that conversions be¬ 
tween observables and metal mass require an assumption of 
physical conditions. Hence it is not straightforward to use 
such inferred measures to discriminate between wind mod¬ 
els. 

In summary, we see that our two wind models distribute 
baryons and metals differently, even though they generally 
give similar observational absorption line properties. This 
suggests that the total amount of CGM gas and metals is 
not sufficiently constrained by the current observations to 
distinguish between competing outflow models. 

4.3 The dynamical state 

A key goal of CGM absorption line work is to be able to 
distinguish whether a given absorber arises from inflowing , 
ou tflowing, or ambi ent material (e.s. iBurchett et al.ll201.T I. 
In lFord et al.l ll2014l b we studied general observational diag¬ 
nostics of the dynamical state of the gas in the ezw model 
and found, for example, that low metal ions tended to trace 
recycled accretion while high metal ions traced ancient out¬ 
flows (both defined below). Here, we expand on this past 
study by examining and contrasting the ezw and cw wind 
models, while focusing on our simulated LOS that mimic 
the COS-Halos data set. In this way we aim to develop in¬ 
tuition about how COS-Halos and similar CGM absorption 
line data sets are able to trace the baryon cycle. 

The first step is to characterise inflowin g and ou tflowing 
gas in the simulations, for which we follow iFbrd et al] (l2014lb 
Briefly, we separate CGM gas particles into five different 
kinematic categories, utilising the past and future informa¬ 
tion we have within the simulations. Since we are interested 
in CGM and not ISM, each of these categories excludes gas 
that is in the ISM at 2 = 0.25. The categories are: 

(i) Pristine Accretion. This is gas that is accreting, 
meaning it is not in the ISM of a galaxy at z = 0.25 but 
will either end up a the galaxy (as an ISM particle or star 
particle) or pass through a galaxy and be ejected as a wind 
particle by 2 = 0. Pristine accretion is accreting gas that 
has not previously been in a wind. 

(ii) Recycled Accretion. This is gas that is accreting, 
as defined above, but has been in a wind prior to 2 = 0.25. 

(iii) Young Outflows. This is gas that was ejected in a 
wind “recently”, relative to 2 = 0.25. We define recently as 
1 Gyr prior to z = 0.25 (z = 0.36 for our cosmology), as this 
is roughly the time a particle would take to leave the halo if 


it were simply launched in a wind, and was not slowed down 
by hydrodynamic interactions. 

(iv) Ancient Outflows. This is gas that was ejected in 
a wind longer than 1 Gyr ago (before 2 = 0.36). 

(v) Ambient. This is gas that will not accrete onto a 
galaxy by 2 = 0, and has never been in a wind by 2 = 0.25. 
One can think of ambient material as gas that is neither 
inflowing nor outflowing, and hence not participating in the 
baryon cycle. 

In Figure [TT1 we plot the EW versus impact parameter, 
similar to Figure [3] except now we decompose the simulated 
equivalent widths into contributions from the five dynamical 
c ategories abo ve. We calculate these in the manner described 
in lFord et al.l (l2014h , by making new simulation outputs con¬ 
taining only those gas particles that fit the definition of the 
various categories, rerunning LOS through them, generating 
spectra, and fitting column densities. The ions are ordered 
from low to high ionisation energies, as in earlier figures. We 
also plot as solid black lines results obtained including parti¬ 
cles from all categories and the ISM. The contributions from 
all the categories (coloured lines) do not necessarily exactly 
add to the total (black line), as they have been fitted sep¬ 
arately. We also show the COS-Halos data for reference, as 
in Figure [3] including detections and limits. To avoid clut¬ 
ter we do not include the observed error bars, which are the 
same as in Figure [3] Finally, we bin both the model and the 
observed points in 25 kpc bins. 

Each of these new simulation snapshots is by definition 
sparser than the complete simulation, so it is no longer ap¬ 
propriate to use the medians of LOS with detections above a 
certain limit, as in Figure[3] To account for the sparser sim¬ 
ulation boxes, in Figure [IT] we now plot the mean of all the 
model LOS at a given impact parameter. As it is not appro¬ 
priate to compare the mean of all LOS (as plotted here) to 
only detections (as in the observations), we caution against 
comparing the solid black line to the observed black points. 
The observed points are shown here to simply guide the eye 
and provide a rough comparison. For a more detailed, fair 
comparison, see Figures [3] [4] and [5] We will now examine the 
behaviour of these ions in the models, with the understand¬ 
ing that these models can have significant discrepancies with 
the observed values, particularly for silicon and Ovi. 

The top panels show HI, in ezw (left) and cw (right). 
The g eneral trend is the same as that identified in lFord et al] 
J2014 ): at low impact parameters, HI predominantly arises 
from recycled accretion, tracing dense gas that is about to 
re-accrete within several Gyr, while at high impact param¬ 
eters it arises mostly in pristine ambient gas. At all impact 
parameters there is a non-trivial contribution from ancient 
outflows, which follows the ambient curve at a level that is 
set predominantly by the ratio of the masses in the ambient 
and ancient outflow categories. These trends are qualita¬ 
tively true in both wind models. 

There are, however, non-trivial differences between 
these wind models. Relative to the ezw model, in the cw 
model the recycled accretion only dominates at the very 
smallest impact parameters. Also at smaller impact param¬ 
eters in the cw model, the contribution from ambient gas 
is 2-3 times more than than that from ancient outflows, 
whereas in the ezw model it is comparable. Despite these 
variations, the total (black) line is very similar in both cases, 
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Figure 11. Equivalent width (A) versus impact parameter for the COS-Halos data set (black symbols both panels), as well as for the 
ezw (left) and cw (right) models, for various ions ordered from low to high ionisation energy. The model points represent the mean value 
of all LOS in 25 kpc impact parameter bins. Broken coloured lines correspond to pristine accretion (PA; dotted light blue), recycled 
accretion (RA; dashed purple), young outflows (YO; dash-dot green), ancient outflows (AO; dash-dot-dot orange) and ambient (AMB; 
long dash red). Solid black lines include all particles that lie along the LOS, regardless of category, and in contrast to Figure [3] they 
include LOS below the detection threshold. 


showing that the HI equivalent width is not able to, on its 
own, discriminate the dynamical state of the CGM gas. 

Our lowest i onisat ion ion is Mgn and, again in accord 
with I Ford et ahl (| 20141 ), we find that recycled accretion is 
the dominant contributor at all impact parameters. Beyond 
100 kpc, there is an increasing fractional contribution from 
ancient outflows, but this is still sub-dominant. This is likely 
because Mgn absorption (and low ion s gene rally) drops off 
very quickly with impact parameter (IFord et al. 120131') : at 
large impact parameters the absorption is likely coming from 
satellites, not the central galaxy. This trend is true in both 
wind models, though cw at small impact parameters shows 
a small contribution from ambient gas and young outflows. 

The trends for Mgn are mimicked by the mid-ion Si IV; 
this ion is still predominantly tracing halo gas that will ac¬ 
crete onto the galaxy. In contrast, O VI shows a dominant 
contribution from ancient outflo ws at all but the smallest 
impact parameters. As argued in IFord et al.1 (120141 '). this is 
consistent with the idea that O VI comes from a diffuse halo 


around the galaxy, enriched over long periods of star for - 
mation (iTumlinson et al] 1201 ll : lOppenheimer et al.l 120121 ') . 
These trends are also fairly similar in the cw versus the ezw 
models, which indicates that these interpretations are ro¬ 
bust at least to the range of wind models that we probe 
here. Note that, as discussed earlier, ours are not represen¬ 
tative of all possible outflow mechanisms in current CGM 
simulations. 

In summary, the most clear trend is that at small impact 
parameters, recycled accretion dominates the absorption of 
low ions (including H i), and that ancient outflows dominates 
Ovi absorption. Most Hi, meanwhile, arises from ambient 
gas at larger imp act parame t ers. These trends mimic the 
expectations from lFord et al] |2014h . Hence the broad intu¬ 
ition developed in our previous work can be used to roughly 
interpret absorption line observations from COS-Halos and 
similar studies, even though they may be less clear owing 
to a conflation of impact parameter and halo mass trends. 
Our two wind models show qualitatively similar behaviour, 


























































































Cycling in Low-Z CGM 19 


indicating that at least within the class of outflow mod¬ 
els considered here these results are robust. The dominance 
of recycled accretion for low ions and ancient outflows for 
Ovi explains why our model predictions for star-forming 
and passive galaxies can be so similar: the absorbing met¬ 
als are predominantly ones that were ejected long ago, and 
they are not coupled to the current star formation rate of 
the central galaxy. Thus, a quenching scheme that shuts off 
star formation but does not directly influence gas already 
in the halo has relatively little impact on CGM observables. 
Similarly, since the ambient gas dominates the HI absorp¬ 
tion, one expects little difference between star forming and 
passive galaxies, except perhaps at small impact parameters 
where recycled accretion becomes important. 


5 SUMMARY 

We have compared our cosmological hydrodynamic simula¬ 
tions including two galactic outflow prescriptions with HI 
and metal absorption line observations in the CGM from 
the COS-Halos survey, to (i) assess how well current models 
that self-consistently enrich the CGM and IGM fare against 
state of the art CGM observations; (ii) determine if available 
observations are able to distinguish between outflow mod¬ 
els; and (iii) provide interpretations about the physical and 
dynamical state of CGM gas as observed by COS-Halos. To 
accomplish this, we extracted LOS around simulated galax¬ 
ies that match the COS-Halos sample closely in stellar mass 
and impact parameter, analysed them in a manner consis¬ 
tent with the observations, and made comparisons to di¬ 
rect observables, namely covering fractions and equivalent 
widths as a function of impact parameter and specific star 
formation rate and the kinematics of the absorbers relative 
to the host galaxy. In particular, we compared to COS-Halos 
data tracing H i, Mg n, Si n, C II, Si iii, Si iv, C ill, and O vi in 
44 sightlines out to an impact parameter of 160 kpc around 
galaxies with stellar masses between 10 9 ' 5 < M» < 10 115 Mq 
at 2 = 0.25. 

Our main conclusions are as follows: 

(i) Our favoured wind model, ezw, is in broad agreement 
with key absorption line observables from COS-Halos for 
Hi and some low ionisation metal lines. In particular, for 
these ions we find good agreement in the observed equiv¬ 
alent widths versus impact parameter, the total equivalent 
width per unit path length, and ion covering fractions above 
the COS-Halos detection limit. This model is also in broad 
agreement with ion ratios and absorber kinematics relative 
to the host galaxy. 

(ii) However, there are numerous discrepancies that hint 
at missing or poorly represented physics in the simulations. 
These include: 

• Both wind models underpredict the two most com¬ 
monly observed ions from COS-Halos, namely the mid-ion 
Sim and high-ion Ovi. For Ovi, the mismatch occurs at 
all impact parameters, while for Si ill the match is in¬ 
creasingly poor towards larger impact parameters. This 
mismatch occurs in both equivalent width per unit red- 
shift and the covering fraction. In general, all the ions 
but O vi have a steeper gradient of absorption versus im¬ 
pact parameter in the models than what is observed, likely 


suggesting that there are not enough metals in the out¬ 
skirts of halos in the models. We also note that both wind 
model simulations have the same resolution, and hence 
fail equally to reproduce any CGM structure below the 
resolution limit. 

• The total amount of absorption in Siiv, a mid ion, 
is generally too high compared to observations, while the 
amount in high ion O vi is too low. The ion ratios Si IV/HI 
and O Vl/H I echo these trends. This is likely to be an issue 
related to how much metals are deposited into different 
phases of the CGM, in the sense that our model puts too 
many metals in ~ 10 4 K gas and not enough in ~ 10 5 K 
gas. This occurs in both wind models, hinting that the 
underlying cause may be the way in which we eject winds 
as cold, unmixing gas from the ISM. 

• The large scatter in the equivalent widths at a given 
impact parameter for low ionisation lines is not well re¬ 
produced in the simulations. This suggests that local ef¬ 
fects perhaps owing to local ionisation sources and self¬ 
shielding might be necessary to fully explain the distribu¬ 
tion of these low ions. 

• While observations show a substantial difference in 
the incidence of OVI absorption around star-forming and 
passive galaxies, our simulations predict similar levels 
of absorption. This discrepancy may arise because our 
phenomenological quenching scheme suppresses star for¬ 
mation in massive galaxies but does not directly affect 
halo gas. Conversely, our models successfully explain the 
observed similarity of low-ion absorption around star¬ 
forming and passive galaxies. 

(iii) The constant wind model produces CGM absorption 
in all ions that is surprisingly similar to that from the ezw 
model, despite substantially different input wind scalings, 
which results in higher halo gas temperatures and fewer cool 
baryons and metals. It is not possible to distinguish the cw 
model from the ezw model using the available equivalent 
width, covering fraction, or kinematic data, although ion 
ratios do reveal some differences. The lack of discrimination 
between the ezw and t he cw models is a l so su rprising given 
our previous results in lFord et al.l ( 2013i . 2014 ) that showed 
larger differences between similar models. We show that this 
owes largely to the fact that previously we compared models 
at a fixed halo mass, whereas here we compare at a fixed 
stellar mass. The CGM metal absorption properties in these 
models are thus more closely correlated with stellar mass 
than halo mass. 

(iv) Both the ezw and cw models show large contribution 
from recycled accretion at low impact parameters (< 50 kpc) 
in low ions including Hi, and from ancient outflows in Ovi 
at most impact paramete rs, in agreement with our previous 
results (|Ford et al.l l2014h These trends are robust for our 
simulations, suggesting that they can be used to roughly 
infer the dynamical state of CGM absorbing gas at least 
within the context of the probed outflow models. 

The agreement in Hi, Mgn, and Cll between hydrody¬ 
namic simulations including galactic outflows and the most 
recent observations of CGM gas offers encouraging support 
for the baryon cycling paradigm. The disagreement between 
both wind models and observations in Si iii and Ovi sug¬ 
gests further work is necessary in the simulations to get 
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outflows to place metals in the proper location and phase 
in the CGM. 

While the COS-Halos observations of equivalent widths 
and covering fractions cannot discriminate between the two 
wind models we consider here, it is important to recognise 
that these span a relatively narrow range of possible out¬ 
flow models. The ezw and cw models have different mass 
loading factors and ejection velocities, but in both cases our 
implementation ejects gas at the ISM temperature (usually 
~ 10 4 K) and does not mix metals from SPH wind particles 
into the surrounding CGM. These assumptions are a limiting 
case, since some mixing must certainly occur, and if super¬ 
novae are involved in driving outflows it is likely that the 
ejected gas will be at least somewhat heated. Other groups 
have utilised fundamentally different algorithms for ejecting 
outflows, for example by super-heating the ISM gas, which 
may give much larger variations in CGM properties than 
either of the models we consider here. We plan to consider 
a wider range of outflow driving mechanisms in the future, 
so this work should be regarded as an initial investigation 
that we hope will set the stage for a more comprehensive 
comparison amongst all the groups working on models of 
the CGM and the baryon cycle. The discrepancies between 
our models and current observations already suggest that we 
must improve our models of outflow ejection, propagation, 
and mixing. 

Within the context of the models we consider here, our 
interpretations of the COS-Halos observations are relatively 
insensitive to the details of how we implement winds. For in¬ 
stance, the conclusion that low ion absorption is dominated 
by enriched gas that will shortly fall into the galaxy, while 
high ions arise in gas that was enriched at earlier times and 
is stably located in hot gaseous halos, is quite generic for our 
models. Of course, these interpretations may be different if 
we fundamentally changed our wind driving mechanisms, 
as discussed above. Nonetheless, the overall consistency of 
our models provides at least a plausible interpretation of 
how to connect the baryon cycle to CGM absorption line 
observations, which will help develop intuition and testable 
predictions for future observational probes of CGM gas. 
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